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In the present work, we investigated the structural and catalytic properties of a prototype system Pt-doped 


CaTiOz by means of first principles calculations. We paid particular attention to the aggregation and 


penetration of Pt on different surfaces of CaTiOz, and subsequent CO oxidation by surface oxygen 


atoms on Pt-doped CaTiOs. Our calculations indicate that CO oxidation can potentially take place when 
Pt is doped on the first layer of CaTiOz(001). The activation barriers are calculated to be 0.20-0.45 eV. 
The possibly induced O vacancy on the surface will produce a magnetic behavior by breaking the spin 


density symmetry due to one Pt-O bond cleavage. Our study is expected to provide an insight into the 


catalytic behavior of Pt ions in Pt-doped perovskite toward the oxidation of exhaust gas. 


| Introduction 


Oxide-supported precious-metal particles have attracted intense 
interest due to their wide applications in heterogeneous catal- 
ysis. However, the well-known strong metal-support interac- 
tions (SMSI) for group VIII (8-10) noble metals on oxide often 
drastically reduce their catalytic performance.’ For example, 
noble-metal-promoted perovskites, potentially used as auto- 
mobile exhaust catalysts, exhibit typical SMSI that leads to 
severe catalyst deactivation.” This deterioration is generally 
compensated by loading an excess of precious metal, which 
could result in overconsumption of limited noble metal 
resources.*” Therefore, effectively maintaining or increasing the 
catalytic activity and controlling the consumption of precious 
metals are highly demanded. Recently, it was observed that the 
incorporation of precious metal into perovskite lattice is helpful 
to increase their catalytic ability. The presence of ionic Pd, Pt or 
Rh in a perovskite matrix, in most cases, is responsible for the 
high activity of noble-metal-promoted perovskites.*® For 
example, BaCeO; was experimentally observed to incorporate 
Pd(u) into B sites to form BaCe,_,Pd,03_5, which maximized 
the catalytic activity for CO oxidation due to the presence of 
cationic Pd(1) in the perovskite host.” Similarly, the experi- 
mentally observed high-temperature activity of LaPdCoO, for 
lean NO, reduction was also linked to the presence of ionic Pd 
in the LaCoO; matrix.’ 

Motivated by the high activity of Pd ions in BaCe,_,Pd,03_5, 
we attempted to investigate the structural and catalytic 
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properties of prototype CaPt,Ti,_,.03 system by means of first 
principles calculations. Actually, the incorporation of noble 
metal into the B site of perovskite (ABO3) lattice has been 
originally found for Pd-doped perovskites (LaFeCoPdO; and 
LaFePdO3) due to the self-regeneration abilities of Pd, i.e., Pd 
could be reversibly incorporated into the perovskite lattice and 
segregate out, as the exhaust gas is cycled between oxidizing 
and reducing conditions.*® The perovskite exhibiting self- 
regeneration feature is named as “intelligent catalysis” (IC), 
which was afterward extended to Pt and Rh-doped CaTiO; 
systems.° By using X-ray and transmission electron microscope 
(TEM) techniques, Scott et al. observed that Rh or Pt in the 
doped CaTiO; could move in and out of the perovskite structure 
with the fluctuation of oxidation and reduction atmospheres. 
In this work, the self-regeneration of single-atom Pt deposi- 
tion on or penetration into CaO- and TiO,-terminated CaTiO; 
system was initially studied to evaluate the possibility of the 
presence of Pt cation. After that, the possible enhancement of 
surface oxygen activity toward CO oxidation by the presence of 
Pt substitution was explored by performing climbing nudged 
elastic band (cNEB) calculations.'°’? The produced oxidation 
product CO, was found to be desorbed from TiO,-terminated 
surface and create an O vacancy, which will induce a spin 
splitting of Pt-5d electron below the Fermi level, to cause 
magnetism transformation. Thus, we expect that the magne- 
tism transition from non-magnet of stoichiometric CaTi,. 
_xPt,O3 to magnetic CaTi,_,Pt,03_; can be applied to monitor 
such oxidation process in the experimental characterization. 


Il Computational methods 


Spin polarized DFT calculations were performed using the 
Vienna ab initio simulation package (VASP). The electron 
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exchange and correlation were treated within the generalized 
gradient approximation (GGA)* using the PW91 functional." 
The projector-augmented-wave (PAW) method of Blöchl* was 
used for the treatment of the core electrons. A soft potential was 
employed for oxygen (Os, 2s’2p’, cutoff energy of 250 eV), while 
a plane wave cutoff of 350 eV was chosen in all calculations. 
Monkhorst Pack (MP) grid’® of (2 x 2 x 1) were used for the 
(2 x 2) CaTiO3 (001) surface. In order to avoid the dipole-dipole 
interactions between periodic slabs, we built the slab models 
consisting of 7 atomic layers with CaO or TiO, as the termina- 
tion for both sides. A vacuum spacing of 14 A was added in the 
perpendicular direction. All atoms were fully relaxed with the 
forces converged to less than 0.03 eV A-+. The Methfessel- 
Paxton technique” with a smearing width of o = 0.1 eV was 
used to minimize the errors in the Hellmann-Feynman forces. 
The optimized lattice constants (a = 5.394 Å, b = 5.512 Å, c = 
7.685 Å) are in good agreement with experimental results (a = 
5.381 Å, b = 5.442 Å, c = 7.641 Å). 

To ascertain whether Pt could penetrate into replace CaTiO} 
host to replace Ti ions, the configurations of CaTi,_,Pt,O3 are 
generated by the cluster expansion (CE) formalism.”*™ For the 
binary system, CE can be written as 

E(@)=Vo+ 5 Visi + 5 Vijojo; + 5 Vinoiojon +... (4) 

i ij ijk 

where V; Vij Vix,..-are the effective cluster interactions (ECIs), ¢ 
is a spin-like occupation variable (e.g., o = 1 if the B-site of 
CaTiO; perovskite matrix is occupied by Ti, and ø = —1 if the 
site is occupied by Pt). The coefficients are obtained by fitting 
eqn (1) to the first principles calculated energies of the selected 
subset from the configurations generated by CE. In principle, 
CE should take all the pairs, triplets, quadruples, etc. into 
account; however, in practice the cluster expansion can effi- 
ciently predict the energy of a system based on fewer coeffi- 
cients. The optimal coefficient can be obtained according to the 
cross-validation (CV) score by using a genetic algorithm.”*** 


Ill Results and discussions 
A. Deposition and migration of Pt on CaTiO;(001) surfaces 


First of all, we attempt to investigate the binding of Pt atom on 
the two commonly observed CaO- and TiO,-terminated 
surfaces. Similar to the previously studied Pd-doped LaFeO; 
system,” the one-fold Oy site at the CaO-terminated surface 
and the four-fold Oy site at TiO,-terminated surface are found 
to be the most favorable binding positions for Pt (see Fig. 1). In 
the calculations, Pt on other trial sites will converge to these two 
sites upon structural optimizations. As shown in Fig. 1a and b, 
Pt binds to the outmost O; on CaO-terminated surface with the 
bond length of 1.93 A and tilting angle (Z PtOO) of about 55°, 
yielding a binding energy (Ep) of 2.93 eV. In the case of TiO,- 
terminated surface, Pt is preferentially located on the Oy hollow 
site, yielding a slightly higher E», of 3.10 eV. Upon the deposition 
of Pt, the four coordinated Oy atoms would be dragged out of 
surface by around 0.30 A to form four Pt-O,; bonds with bond 
lengths in the range of 2.01-2.05 A. 

Interestingly, the charge distribution analysis displays some 
distinctive features for the charge transfer between the two 
surfaces. Specifically, Pt gains 0.55 electron from the CaO- 
terminated surface, whereas donates 0.63 electron to the TiO,- 
terminated surface. This can be understandable since the CaO- 
terminated and TiO,-terminated surfaces are electron-rich and 
electron-poor, respectively. According to the Bader charge 
analysis, the net charges on the CaO-terminated and TiO,- 
terminated surfaces are calculated to be —0.30|e| and +0.73]e|, 
respectively. In the Bader analysis, each atom is divided based 
on the zero flux surfaces those are 2-D surfaces on which the 
charge density is a minimum perpendicular to surface. To 
present further evidence of electron transfer between Pt and 
CaTiO; surfaces, we also calculated the electron density differ- 
ence between the adstructures and separate components with 
each atom frozen at their optimized positions. It is clearly 
shown in Fig. 1c and c’ that Pt-O bonds are formed on both 


Fig. 1 Single Pt atom adsorption on CaO and TiO2-terminated surface. (a-d) provide the topview and sideview of single Pt atom deposition on 
CaO-terminated surface, and its 2D and 3D electron energy difference; (a’—d’) shows the corresponding parts of single Pt atom deposition on 
TiOz-terminated surface. The green, grey, red and blue balls represent Ca, Ti, O and Pt atoms, respectively. 


surfaces after electron transfer. The 3D yellow spherical shapes 
(denoting electron-rich area) in Fig. 1d indicates that Pt serves 
as electron acceptor on the CaO-terminated surface. In contrast, 
the light-blue shapes (denoting electron-poor area) of Pt on 
TiO,-termiated surface shows electron donating behavior. 
Since metal aggregation is usually the major concern for the 
supported catalysts, we next studied the aggregation tendency 
of Pt on the two surfaces. As only one favorable binding site is 
identified for each surface, we evaluated the migration of Pt 
between two nearest Or sites on CaO-termination and two 
nearest Oy, sites on TiO,-termination. The transition state is 
identified by performing cNEB calculations along the 
prescribed pathways, as shown in Fig. 2. The calculated energy 
barriers for Pt migration on CaO and TiO,-terminations are 1.18 
and 1.71 eV, respectively. Such high barriers indicate that the 
migration of Pt is difficult to take place between two nearest Or 
or Oy sites. An intermediate state of O-Pt-O bonding is found 
on the migration pathway of Pt/TiO,-term, which is also 
obtained by the initial optimization as a meta-stable state. 


B. The incorporation of Pt into CaTiO, lattice 


We next estimated the possibility of Pt penetration into the B 
site of CaTiO; lattice in terms of thermodynamics. Here, we 
calculated the binding energy (Ep = Etotat — Evacancy-ti — Ept) for 
Pt with the Ti vacancy at different layers, as shown in Fig. 3a. Ep 
is calculated to be 9.52 eV for the substitution of the 1*-layer Ti 
(denoted as "Pt;;'*" hereafter), which is increased to 9.81 eV 
when Pt replaces the second (denoted as "Pt;;7"“" hereafter) and 
third layer Ti (denoted as "Pt;;°""" hereafter). The larger Ep of 
Ptr”? and Pt;;"4 means Pt tend to penetrate deep into the 
CaTiO; lattice, which is consistent with the experimental 
observation that Pt could dissolve into the CaTiO; lattice under 
oxidizing atmosphere.® Actually, it is reasonable that E» for 
Pt;;"" is greater than that of Ptr: because Pt forms five Pt-O 
bonds in a square-pyramid coordination with the bond length 
of 1.97-2.03 A in Pty“ while it forms six Pt-O bonds in octa- 
hedron coordination with the bond length of 2.02-2.03 A in 
Pty nt. 
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Fig.2 The migration barriers of single Pt atom between the nearest O; 
and Oy sites on CaO- and TiO2-terminated surface, respectively. 
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Fig. 3 (a) The structures of Pt penetration into different layer Ti to 
form Pt-—CaTiOz; (b) convex hulls of the formation energies per 
formula unit (FE/f.u.) of Pt doped CaTiOz structures based on first 
principles calculation and cluster expansion. Data obtained by VASP 
and CE are colored red and green, respectively. 


To further ascertain whether Pt could replace Ti ions of 
CaTiO; and the corresponding stability, we calculated the 
formation energies of CaTi,_,Pt,O3; by DFT methods according 
to Er = Ecrpo — XEcpo — (1 — x) Ecro, where Ecrpo,; Ecpo and 
Ecro are the total energies of CaTi,_,.Pt,03, CaPtO3 and CaTiO3;, 
respectively, x is the Pt doping concentration. The formation 
energies obtained by DFT calculations were verified by the CE 
prediction and a good correlation between them was obtained. 
The weighted CV score is 3.3 meV/f.u. with a root mean square 
(rms) error of 1.4 meV/f.u. in the CE validation process. The 
comparison between the DFT calculated (red points) and CE 
predicted (green points) formation energies are shown in 
Fig. 3b, respectively. The convex hull obtained by CE (green line) 
is well-converged and accurately reproduces that of the DFT 
ground states (red line). The negative formation energies of 
CaTi,_,Pt,O3 indicate the feasibility of Pt occupying the Ti sites. 
However, the energy differences between the ground states and 


other structures are only several meV/f.u., implying that the 
CaTi,_,Pt,O3 structures would probably become a disordered 
solid solution above moderate temperatures. 


C. CO oxidation by surface oxygen 


Since the lattice oxygen has been verified to not oxidize CO in 
bare TiO,(110) surface,** it is natural to ask whether such 
oxidation process will occur in the presence of ionic Pt in TiO,- 


terminated surface. To explore the possible catalytic 
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enhancement by Pt incorporation into the perovskite host, we 
focus on the CO oxidation by surface oxygen in Ptr% and 
Ptr ®t TiO,-termination. Inspired by the recent observation of 
single Pt atom loaded on FeO, nanocrystallite support or gra- 
phene to form Pt,/FeO,, or Pt/graphene, respectively,*”** we built 
a model of single Pt atom deposited on the surface or occupying 
the B site of CaTiO,; to form CaPt,Ti,_,03; to minimize the 
computational cost and examine its catalytic behavior. We 
found that CO can be favorably adsorbed on the Pt site to form 
linear chemisorption state via the C or O ends in the Pty; t 
structure (see Fig. 4a). Apparently, the C end can interact with 
the ionic Pt more strongly, with an adsorption energy (Eads. = 
Ecreo + Eco — Eco-crro) of 1.59 eV and shorter Pt-C distance of 
1.86 A. In contrast, the O end can interact with Pt with an Faas of 
0.16 eV and O-Pt distance of 2.13 A, respectively. For compar- 
ison, we also investigated the CO adsorption on Pt deposited 
onto stoichiometric TiO,-terminated surface instead of replac- 
ing the Ti ion. Upon optimization, the CO molecule initially 
placed 2.1 A above the surface Pt atom, would be pushed away 
in spite of C and O ends, as shown in Fig. 4c, which illustrates 
that the Pt atom deposited on TiO, surface completely loses its 
reactivity toward CO molecule. Actually, this is reasonable 
because the four lobes of Pt-5d orbital become inaccessible after 
the formation of four Pt-O coordination on TiO, surface, 
similar to our previous study of H, chemisorption on Pt/MoO3;.** 

For the Ptr”? structure, the Pt atom is embraced under- 
neath and not available for interacting with CO. Nevertheless, 
CO can be still linearly chemisorbed on the Ti site via the C end. 
The distance between C end and Ti is calculated to be 2.34 A, as 
depicted in Fig. 4b. Although the corresponding Eaas = 0.10 eV 
is significantly weaker than that of Ptr; *, the reactive activity of 
Ti is actually enhanced by loading Pt since CO can be hardly 
adsorbed on the TiO,-terminated surface of the pure CaTiO3. 


pees 


\b-— #-— M 


Fig. 4 The bottom presents adstructures of CO adsorption on Ptr, 


Pty 2"? and Pt deposition on TiO2-terminated surface through C end. 
The top is the corresponding calculated 3D electron density difference 
of Pt-doped CaTiOz(001) with TiO2-termination. 


The similar phenomenon can be found in the case of CO acti- 
vation on BaCeo.soPdo.1003—-5 and undoped BaCeOs, in which 
the higher activity associated with BaCeg 99Pdo.1903_5 is reflec- 
ted by lower light-off temperature (T5o = 83 °C), compared with 
Tso = 225 °C for the undoped BaCeO3.” 

To ascertain the intrinsic mechanism of the improved 
surface activity, the 3D electron density difference is plotted for 
these Pt-doped and undoped structures in Fig. 4. For the Pt;;'* 
structure, it is clearly shown that the five-coordinated Pt atom 
accumulates some electrons (blue contour) along the surface 
normal (Fig. 4a) due to the absence of one Pt-O bond. The 
remained electron on the unsaturated Pt atom is responsible for 
the high reactivity toward CO adsorption. After further incor- 
porating Pt into the second TiO, layer, the surface Ti atom 
obtains some electrons through the Pt-O-Ti bond (Fig. 4b), 
which enhances the reactivity of Ti compared to other surface Ti 
ions. In contrast, no electron accumulation along the surface 
normal can be observed nearby the Pt atom deposited on the 
stoichiometric TiO,-terminated surface. As shown in Fig. 4c, 
only yellow contour of denoting electron deficient in this 
direction was remained due to the tetra- Pt-O bonds formation. 

As discussed above, TiO,-terminated surface exhibits elec- 
tron-poor feature evidenced by Bader charge analysis, which 
results in the ability of attracting electron in the surface. We 
next evaluate the possibility of CO oxidation on surface oxygen 
following the two chemisorption states. For the adstructure 
(Pt-CO) through C end, the C atom is able to further bind with 
the neighbored O, and drags it out of the surface by about 0.90 
A to form a bent CO, structure, yielding a bent angle Z OCO = 
140° and E,as of 2.08 eV. We term this state as CO,-bending. 
From the kinetic point of view, the adsorbed CO needs to 
overcome a barrier of 0.45 eV at the transition state TS; to bind 
with O, atom, as depicted in path I of Fig. 5. The bond length of 
two newly formed C-O bonds is slightly elongated from 1.18 A 
to 1.21 A due to the slight disturbance of C=O triplet bond. 

The oxidation process is also found to potentially take place 
for CO adsorbed through the O end (Pt-OC). Since adsorption 
on Pt’ via the O end results in smaller energy, the upstanding 
C atom will continue to increase its E,g, to 0.34 eV by binding 
with O, atom with a bent angle of Z OCO = 106°. In contrast to 
the n°(CO) structure, this bent CO,-TM state should be ascribed 
to n°(COO) state because the two O atoms are directly bonded 
with Pt. The newly formed C-O, bond of 1.38 A will also drag the 
O, atom out of TiO, surface by about 0.8 A. The calculated 
activation barrier of 0.20 eV toward CO oxidation is even smaller 
than that of the commonly used Au-Cu catalyst (activation 
barrier = 0.60 eV).*° In the TS,', the original C=O bond length 
is elongated to 1.20 A and C-O, distance is decreased to 2.34 A. 
More interestingly, such oxidation state n*(COO) is actually an 
intermediate and would converge to the most stable CO,- 
bending state 7?(CO) eventually. The similar conversion has 
also been found in the previous CO,-TM state, in which COO 
binding is not a stable state.***’ In the TS," (Fig. 5), the distance 
of C-Pt is shortened to 2.67 A, while the original Pt-O bond is 
increased to 2.48 A, yielding an activation barrier of 0.26 eV. 

Whether the product (free CO, molecule) could be desorbed 
from such 7’(CO)-Pt complex is another key issue for the 
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Fig. 5 The CO oxidizing process by the surface O1 atom on TiOz- 
terminated surface of Pt-substitution the first layer Ti. Path | and II refer 
to the two reaction pathways through C and O ends. The corre- 
sponding structures and the bond length are also presented. 


catalytic activity. As depicted in Fig. 5, the CO, desorption will 
increase its thermal energy by about 1.45 eV compared to 
n’(CO)-Pt complex, yielding an activation barrier of 1.45 eV in 
this process. Such a high activation barrier is generally difficult 
to conquer, however, the exothermic energy (about 2.1 eV) of 
forming 7*(CO) in path II possibly accumulate enough energy to 
enable the desorption of CO,. It is noteworthy here that the 
high-barrier reaction pathway II only plays a minor role in the 
CO oxidation process, compared to the low-barrier reaction 
pathway I. Furthermore, an O-vacancy on TiO,-terminated 
surface will be induced by CO, desorption, as shown in the 
dashed box in Fig. 5. Although no obvious deformation can be 
discerned except for the slightly elongated bond length of Pt-O 
bond opposite to Pt-Oyacancy, the electron feature near Fermi 
level is significantly different from the stoichiometric Pt- 
CaTiO . Fig. 6a clearly shows that the whole electronic bands 
are downshifted by approximately 0.5 eV in TiO,-terminated 
surface with an O-vacancy. An obvious spin splitting at valence 
band maximum (VBM) is discerned by destroying the 
symmetric electronic bands of stoichiometric CaTi,_,Pt,O3, 
although the band gap is nearly unchanged. Such electron 
variation will probably induce magnetism transformation from 
non-magnet of stoichiometric CaTi,_,Pt,O3 to magnetic CaTi,. 
_xPt,O3_5, which is further confirmed by the calculated total 
magnetic moment of 1.90 ug in CaTi,_,Pt,O3_5. Projected 
density of states (PDOS) indicates that the split electronic band 
in VBM is mainly ascribed to the splitting of Pt-5d bands, as 
shown in Fig. 6b, which is consistent with the calculated local 
magnetic moment of 1.10 ug on Pt. Compared to the original 
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Fig. 6 Comparison of the total density of states (TDOS) (a) for O- 
vacant (Pt—CaTiOz_,) and Pt-CaTiOz and the projected local density of 
states (PDOS) of Pt-5d band (b) to illustrate their electromagnetic 
behavior near Fermi level. 


smaller peak, Pt-5d band becomes more narrow and stronger. 
Such emergence of magnetism is contributed by the broken 
symmetric square-pyramid structure of five-coordination of Pt, 
and is expected to be helpful to monitor the surface oxygen 
function states in the CO oxidation process. 


IV Conclusion 


Based on the DFT and CE calculations, both the increased 
binding energy with Pt penetration and the negative formation 
energy of CaTi,_,Pt,O3 indicate the feasibility of replacing Ti 
sites in CaTiO3 by Pt, which is consistent with the experimen- 
tally observed self-regeneration of Pt particles on CaTiO3. The 
penetrated Pt ion could enhance the surface activity of CaTiO; 
to form two adsorption states through the C and O ends. 
Further, we show that following oxidation by the surface O atom 
could readily take place to yield more stable n?(CO)-Pt complex, 
due to the low activation barrier of 0.45 eV and 0.20 eV in the 
pathways I and II. From the whole reaction energy profile, CO 


oxidation path I is dominate process owning to the lower 
thermal energy along the oxidation process. However, the minor 
oxidation path II is possible to induce CO, desorption and 
hence cause a magnetic transition (ie., from non-magnetic 
CaTiO; to magnetic CaTiO3_ 5). Our study is expected to provide 
an insight of catalytic behavior of Pt ion in Pt-doped perovskite 
toward exhaust gas oxidation. The magnetism transition could 
possibly used as a measurement to monitor the oxidation 
process. 
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